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A  METHOD  FOR  CALCULATING  THREE-DIMENSIONAL  FLOW  IN  BLADED  TURBOMACHINES  BY  USING 
ARBITRARY  QUASI  PERPENDICULAR  INTERSECTING  PLANES 

Xin  Xiaokang,  Jiang  Jinliang 

Fluid  Dynamics  Section,  Department  of  Mathematics,  Fudan  University 

Since  the  Great  Proletariat  Cultural  Revolution,  considerable  development 
has  been  seen  in  the  research,  design  and  manufacture  of  bladed  turbomachines  in 
China.  To  further  increase  the  efficiency  and  power  output  with  small  dimensions, 
at  present  the  theory  of  three-dimensional  fluid  flow  is  being  applied  to  upgrade 
design  calculations  in  many  plants  and  organizations  of  scientific  research.  We 
follow  Great  Leader  Chairman  Mao's  instruction,  "Revolution  in  Education"  by 
meeting  the  demands  of  industrial  and  agricultural  production  as  a  research 
subject.  Along  with  the  Shanghai  Shipyard,  Shanghai  Diesel  Engine  Plant,  and 
Shanghai  Institute  of  Internal  Combustion  Engines,  we  are  engaged  in  the  design 
and  calculation  of  three-dimensional  fluid  flow  in  superchargers.  This  paper 
presents  one  result  in  our  open-door  educational  activity. 

Reports  (1,2)  first  proposed  fundamental  equations  and  calculation  processes 
in  the  three-dimensional  fluid  flow  in  bladed  turbomachines.  This  calculation 
method  aims  at  simplifying  the  three-dimensional  flow  field  in  a  bladed  turbo¬ 
machine  into  two  two-dimensional  flow  fields  (S^  and  flow  surfaces);  by 
mutually  successive  substitutions  of  these  two  two-dimensional  flow  fields. 

Translator's  note;  flow  surface  or  surface  of  flow  should  read 
throughflow  surface. 
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solutions  of  thrce-d imons iona 1  flow  field  can  be  obtained.  Even  at  the  present 
time,  this  calculation  method  is  widely  adopted  both  domestically  and  abroad. 
However,  because  a  great  volume  of  calculations  is  required  to  mutually  and 
successively  substitute  between  the  two  groups  of  flow  surfaces,  the  present 
calculation  method  (3-10)  of  three-dimensional  fluid  flow  does  not  employ  succes¬ 
sive  substitutions.  Instead,  one  successive  substitution  from  the  to  flow 
surface  is  used  to  derive  the  quasi  three-dimensional  solution. 


From  the  method  of  arbitrary  quasi  perpendicular  intersecting  lines  (6-8)  in 
the  streamline  curvature  method,  we  developed  a  method  of  arbitrary  quasi 
perpendicular  intersecting  planes.  In  calculations,  the  mutual  influence  between 
the  two  groups  of  flow  surfaces  is  simultaneously  considered,  so  that  no  corrections 
need  to  be  made  from  each  other  by  successive  substitution  between  these  two 
groups  of  flow  surfaces.  In  this  method,  the  theory  of  a  variable  curve  sample 
is  applied.  After  appropriate  successive  substitutions  of  spatial  streamlines, 
the  three-dimensional  solution  of  the  entire  blade  channel  can  be  directly  calcu¬ 
lated,  including  shapes  of  arbitrary  and  S.,  flow  surfaces,  shapes  of  spatial 
streamlines,  and  the  distribution  of  flow  parameters,  such  as  velocity  and  pressure. 

As  the  first  step  in  solving  the  real  three-dimensional  flow  field,  we  apply 
the  approximate  hypothesis  that  the  flow  surface  is  a  surface  of  revolution; 
this  hypothesis  is  applied  in  most  reports.  Thus,  there  are  briefer  and  more 
convenient  steps  in  flow  integration  and  reverse  interpolation  of  equal -sector 
flow  into  new  streamlines.  Hence, the  solution  obtained  is  basically  still  the 
quasi  three-dimensional  solution.  The  next  step  is  to  eliminate  the  hypothesis 
that  the  surface  is  a  surface  of  revolution  in  order  to  obtain  the  solution  of 
a  real  three-dimensional  flow  field.  Then  the  curvature  of  flow  surface  can 
be  calculated. 

I.  The  Method  of  Arbitrary  Quasi  Perpendicular  Intersecting  Planes 

In  principle,  this  method  can  be  applied  to  axial-flow,  diametrical -flow, 
and  mixed-flow  bladed  turbomachines.  The  paper  uses  an  example  of  a  bladed 
channel  in  a  centrifugal  compressor  bladed  turbomachine  for  a  detailed  illustration. 

Figure  la  shows  a  bladed  channel  of  a  centrifugal  compressor  bladed  rotor. 
Figure  lb  shows  the  shape  of  the  meridian  plane  of  the  bladed  channel.  The 
report  (6)  introduces  the  situation  of  utilizing  an  arbitrary  quasi  perpendicular 
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interesting  line  qm  to  solve  for  meridian  flow  field  (S,  flow  surface).  The  report 

(7)  introduces  a  situation  in  which  use  is  made  of  an  arbitrary  quasi  perpendicular 

intersecting  line  q^  to  solve  for  the  trans-blade  flow  field  (S^  flow  surface). 

In  the  ordinary  method  of  the  quasi  perpendicular  intersecting  line,  first  we 

integrate  the  velocity  gradient  equation  in  the  direction  of  q  or  q,..  Next,  bv 

m  U 

flow  correction  the  velocity  distribution  along  the  quasi  perpendicular  intersecting 
line  can  be  obtained.  Then,  by  inverse  interpolation  of  equal-sector  flows,  the 
position  of  new  streamlines  at  the  S£  or  flow  surface  can  be  obtained.  In  the 
paper,  the  quasi  perpendicular  intersecting  lines  can  constitute  a  curved  surface 
Sq  (See  Fig.  la),  which  is  formed  by  the  quasi  perpendicular  intereseting  straight 
line  qm  (Fig.  lb)  revolving  around  the  z  axis  for  a  complete  circle.  We  then  take 
this  curved  surface  as  the  quasi  perpendicular  intersecting  plane  in  solving  the 
problem.  First,  the  fundamental  equations  can  be  used  to  derive  a  velocity 
gradient  equation  along  an  arbitrary  direction  in  space.  Then  two  velocity 
gradient  equations  can  be  derived  along  the  q  and  q„  directions.  Integrate  these 
two  equations.  Next,  through  flow  correction,  the  velocity  distribution  at  the 
quasi  perpendicular  intersecting  plane  can  be  obtained.  Then,  by  reverse  interpola¬ 
tion  of  equal-sector  flow,  new  positions  of  spatial  streamlines  can  be  obtained. 


a»(6) 


Fig.  1. 

Key:  (1)  Pressure  surface;  (2)  Quasi  perpendicular  intersecting 
plane;  (3)  Rotor;  (4)  Stator;  (5)  Suction  surface;  (6)  Streamline; 
(7)  Quasi  perpendicular  intersecting  line. 

The  fundamental  hypothesis  adopted  by  the  paper  is  as  follows: 


1.  The  working  medium  (or  fluid)  is  a  perfect  gas,  which  is  viscousless  with 
equal  C^.  The  flow  is  isentropic.  Friction  losses  in  a  real  gas  can  be  adjusted 
by  total  pressure  loss  as  indicated  in  reports  (6,7),  or  revealed  by  using  the 


ordinary  polytropic  process. 


2.  The  flow  passing  through  the  blade  rotor  is  of  finite  length,  relative  to 
blade  rotor. 

3.  The  Sj  flow  surface  is  a  surface  of  revolution.  This  hypothesis  is  iden¬ 
tical  to  the  hypothesis  used  in  most  reports,  and  restricts  the  application  range 
of  this  method. 

We  now  proceed  to  derive  the  fundamental  equations  used  in  the  paper. 


In  reports  (6,7),  the  velocity  gradient  equation  at  any  curve  direction  in 
the  space  along  the  flow  field  is  as  follows: 


rfir  d,  ,  .  dx  ,  dd 
—  a -  +  b  —  -T  e  — — 


In  equation  (1), 
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Here,  W  is  the  relative  velocity;  q  is  the  intercept  of  the  space  curve;  r£  is  the 
radius  of  curvature  of  the  projection  by  the  streamline  on  the  meridian  plane 
(briefly  called  the  meridian  streamline);  a  is  the  included  angle  between  the  z  axis 
(rotating  axis)  and  the  tangent  of  the  meridian  streamline;  (I  is  the  included  angle 
(Fig.  2a, b)  between  the  meridian  plane  and  the  relative  velocity  vector;  \  is  the 
previous  rotation  (A=rVy);  h!  is  the  total  entropy  of  each  streamline;  and  is 
the  angular  velocity  of  the  blade  rotor  rotation. 

From  Fig.  2,  the  following  geometric  relationship  equation  can  be  established: 
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Here,  Wfl,  and  are,  respectively,  the  relative  velocity  components  along 

the  diametrical  direction,  circumferential  direction,  and  axial  direction:  W  is 

in 

the  relative  velocity  component  of  meridian  plane;  dm  is  the  micro-element  of 
the  arc  length  of  the  meridian  streamline. 


Fig.  2 

Key:  (1)  Rotor;  (2)  Stator;  (3)  Streamline;  (4)  Quasi 
perpendicular  intersecting  line. 


In  the  appendix,  velocity  gradient  equations  can  be  derived  in  the  direction 
of  two  groups  of  curves  along  the  quasi  perpendicular  intersecting  plane: 


One  velocity  gradient  equation  is  in  the  direction  of  the  intersecting  line 
between  the  quasi  perpendicular  intersecting  plane  and  the  blade  pressure  surface 
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Here,  <//  is  the  included  angle  between  the  z  axis  and  the  normal  line  n^  of  the 
intersecting  line  AB  referred  to  concisely  as  the  quasi  perpendicular  intersecting 
line  between  the  meridian  plane  and  the  quasi  perpendicular  intersecting  plane; 
qm  is  a  parameter,  indicating  the  intercept  along  line  AB  starting  from  the  inter¬ 
secting  point  c  between  line  AB  and  the  rotor;  0<j=f<j(r,z)  showing  the  pressure 
surface  equation  of  the  blade. 

Another  velocity  gradient  equation  is  in  the  arc  direction  (0  direction)  when 
qm=constant  along  the  quasi  perpendicular  intersecting  plane. 

-  CU' +  O  +  —  (6) 

je  wfl  dd  /  1  ’ 

Here,  C  and  D  are  the  same  as  in  Equation  (5). 

The  continuous  equation  is  satisfied  by  using  the  same  flow  (equivalent  to 

a  given  flow)  passing  through  the  quasi  perpendicular  intersecting  plane.  The 

following  is  the  equation  of  flow  weight  as  the  flow  passes  through  an  arbitrary 

quasi  perpendicular  intersecting  plane  and  stays  within  the  channel  of  two  adjacent 

blades.  r«»,  ... 

G  •*  ]  pg W  ■  n  •  rdOdqm 

—  I,"*  —  a)rd8d<]m  (7) 

i'.i  Equation  (7),  0^  is  the  angle  of  the  blade  pressure  surface  and  0g  is  the  angle 
of  the  blade  suction  surface.  As  integrations  are  conducted  along  a  curved  surface, 
0d*Hj(qm)  and  8s=0s(qm);  these  are  equations  of  the  curve. 

These  three  equations  are  the  fundamental  equations  we  adopted. 

II.  The  Solution  Procedures  in  This  Method 

The  procedures  in  this  method  are  as  follows: 


1.  First  on  the  meridian  plane,  arbitrarily  select  m  quasi  perpendicular 
intersecting  lines  (Fig.  lb),  which  are  quasi  perpendicularly  intersecting  with 
all  meridian  streamlines.  For  sake  of  simplicity,  generally  straight  lines  are 
selected.  Then  the  corresponding  angle  0  of  these  m  quasi  perpendicular  inter¬ 
secting  lines  can  be  written  as  /  *  —  \ 


0,  «  *rc  »in 


/  *U  -  «» 

*uY  +  r,y 


0  -1,  2, 


.)  (8) 


6 


In  the  equation,  (r^,  z^)  and  (r.,.,  are,  respectively,  the  coordinates  of 
the  intersecting  points  between  the  quasi  perpendicular  intersecting  line  on  one 
side,  and  the  rotor  and  stator,  on  the  ether. 

2.  Rotate  this  group  of  quasi  perpendicular  intersecting  lines  about  the  z 
axis,  we  obtain  the  required  quasi  perpendicular  intersecting  planes. 


3.  Form  the  initial  meridian  streamlines  on  a  meridian  plane.  As  usual,  we 
can  adopt  an  equidistant  division  of  quasi  perpendicular  intersecting  straight 
lines,  or  make  the  revolving-ring-shaped  channel  areas  equal  in  value.  We  adopt 
the  first  method  by  dividing  into  n  equal  parts  to  obtain  (n+1)  meridian  stream¬ 
lines.  Then  we  can  see  that  the  coordinates  (rD,  z  )  of  the  intersecting  point 

D  D 

between  the  initial  meridian  streamlines  and  the  quasi  perpendicular  intersecting 
lines  are :  .  _ 

“  r„  +  - (r*  —  r„) 
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4.  In  the  hypothesis,  the  flow  surface  is  a  surface  of  revolution, 
therefore,  as  mentioned  above,  revolving  the  (n+1)  meridian  streamlines  around 
z  axis  will  generate  (n+1)  flow  surfaces.  The  initial  space  streamlines  can  be 
generated  on  (n+1)  flow  surfaces. 

First,  the  cutting  of  a  blade  channel  of  the  rotor  by  a  quasi  perpendicular 
intersecting  plane  is  drawn  in  Figure  3.  Figure  3a  is  the  developed  diagram 
(projection  on  a  plane  with  z=constant)  in  an  ordinary  rectangular  coordinate 
system.  Figure  3b  is  the  developed  diagram  in  the  qm-0  coordinate  system. 

From  Figure  3‘,  when  a  quasi  perpendicular  intersecting  plane  cuts  through 
the  flow  channel,  the  boundary  of  the  cross  section  thus  produced  forms  a  curved 
parallelogram.  On  qm-0  diagram,  the  upper  and  lower  sides  of  the  parallelogram 
are  straight  lines:  V°  and  Vqm  stator 

Moreover,  the  intersecting  line  between  the  flow  surface  and  the  quasi  perpen¬ 
dicular  plane  is  also  a  straight  line: 

1~  ^ ’^H^m  stator 
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Therefore,  along  these  straight  lines,  r=constant  and  :=constant.  The  var  ation 
range  of  angle  0  is  64<e<£6, 


(•)  (b> 

Fig-  3 

Key:  (1)  Pressure  surface;  (2)  Suction  surface. 

Now,  we  make  h  equal  divisions  on  0  at  the  straight  lines  qm=constant ;  then 
we  have  (n+l)(h+l)  initial  space  streamlines.  Then  we  can  derive  the  result  that 
the  coordinates  of  0  of  the  intersecting  point  between  the  quasi  perpendicular 
intersecting  plane  and  the  initial  space  streamlines  are: 

(»'  *“  1 .  2 ,  •  •  • ,  m  \ 

1,  2,  •••,  »  +  1  (10) 
K  -  1.  ?,•••,  h  +  1/ 

In  the  equation,  ©^(r.z)  and  ©s(r,z)  are  determined  by  the  blade  forming  equation 

of  air-guiding  rotor  and  working  rotor.  The  analytical  equation  can  be  derived 

or  values  of  divergent  points  can  be  obtained.  Thus,  after  r  and  z  are  given, 
values  of  ©^  and  qs  can  be  solved  or  obtained  by  interpolation. 

Thus,  we  obtain  nodes  of  initial  streamlines  in  the  entire  blade  channel. 
There  are  altogether  m  quasi  perpendicular  intersecting  planes,  (n+1)  flow 
surfaces  and  (h+1)  S2  flow  surfaces.  Coordinates  of  intersecting  points  (nodes) 
of  these  three  groups  of  surfaces  can  refer  to  Equations  (9)  and  (10). 

5.  After  obtaining  the  coordinates  of  the  initial  streamlines,  according  to 

the  theory  of  sampling  curve  we  can  obtain  various  geometric  quantities  (such  as 
«,  (1,  and  rc)  of  streamlines  and  solve  for  various  parameters  (such  as  dh'^/dm 
and  dW^/dm)  required  in  Equations  (4),  (6)  and  (7).  The  solution  procedures  are 

8 


the  sane  as  in  report  (A) .  As  in  report  (6),  the  initial  value  of  W  can  be 
assumed  constant;  in  our  case,  the  initial  value  of  W  arc  assumed  as  constant  at 
each  quasi  perpendicular  intersecting  plane  for  more  rapid  convergence  in  succes¬ 
sive  substitutions.  Here,  the  initial  value  of  is  given  by  input. 

6.  Same  equations  for  numerical  solution  can  be  adopted  in  the  numerical 
integration  of  two  differential  equations  (4)  and  (6).  We  still  adopt  the 
prediction  -  correction  method  (i.e.,  Runge-Kutta  method)  described  in  report  (6). 

In  this  method,  the  procedure  for  solving  W  values  of  the  entire  flow 
field  is  as  follows:  To  each  quasi  perpendicular  intersecting  plane,  first  solve 
for  W  values  (Fig.  3b)  at  points  A^ ,  A^,...,  A^+^  with  initial  values  of  point  A^ 
by  using  pressure  surface  integration  equation  (4).  Then,  use  qm=constant  to 
integrate  linear  differential  equation  (6),  from  point  A^  to  point  at  the 
suction  surface  in  order  to  obtain  W  values  of  the  entire  net  lattice  points. 


7.  Use  flow  conservation  equation  (7)  to  check  whether  the  total  flow  can 
be  satisfied  in  order  to  revise  the  initial  value  of  W  for  point  A^ . 

CK  —  (  ’  ptW  COSTCO s(</>  —  a)J6  (11) 

J*4 

Then  the  flow  weight  equation  (7)  passing  through  any  arbitrary  quasi  perpendicular 
intersecting  plane  can  be  written  as 


(*»  stator 


*  (r,  +  qmcOttJ>)Jqm 


These  integrations  can  be  easily  obtained  by  using  the  sampling  curve  theory  of 
the  divergent  point.  Moreover,  a  group  of  successively  additive  numbers  can  be 
obtained.  If  the  calculated  total  flow  G^  and  the  given  flow  Gq  do  not  satisfy 
precision  requirement  \G,  —  Gt\  <  ec  (,cc  input)  (13) 

Then,  we  can  use 

-WA,&  (14) 

G, 

to  be  the  new  initial  value  of  the  velocity  at  point  A^.  Repeat  steps  6  and  7 
until  the  precision  requirement  (13)  is  satisfied.  Generally,  this  process 
requires  five  to  six  iterations. 


8.  After  satisfying  the  precision  requirement  (13),  apply  the  method  of 

reverse  interpolation  of  equal -division  flow  and  solve  for  the  streamline 

y 


coord i nat  os  (  r,  , 

Vi 


U 


and  o.l. 

v» 


First,  apply  reverse  interpolation  by  n  equal  divisions  of  meridian  flow  at 

the  total  flow  curve  G.-q  (fig.  4a)  and  solve  for  various  values  of  (q  )  ,  then 

i  in  'm  G 


rc  ~  r,  +  ($«,)<;  •  cos  4> 
*c  “  *i  —  (?«)c  •  sin  <l> 


(15) 


Apply  reverse  interpolation  to  angle  0  and  reverse  interpolation  for  h  equal 

divisions  at  the  flow  curve  G,-0  (Fig.  4b),  we  obtain  the  value  of  o  .  However, 

K  G 

since  values  of  r  and  z  have  changed,  values  of  and  0^  for  boundary  points 

should  start  from  r.  and  z^  values  using  blade  equations  0^=0(j(r,z)  and  0  =9s(r,z) 

a:. 4  solving  for  values  of  (6  )  and  (9  )  .  Here  we  employed  an 

a.  G  S  G 

approximation  hypothesis:  we  considered  that  the  9  values  obtained 

G 

for  the  original  r  and  z  can  be  approximately  regarded  as  the 

partition  values  for  r  and  z „  but  0O  and  9  nonetheless  use  the 

G  G  cl  s 

values  (9^)^  and  (0s^g'  Because  the  loosening  factor  n  is  adopted 
later  on,  this  approximation  will  not  affect  convergence. 


Fig.  ^  . 


9.  Lastly,  we  derive  the  maximum  deviation  between  new  and  old  streamline 

Pil  I  nt  S 

A/—  max  K'c  *”  riY  +  (*c  ““  *»)J 
!<*<"• 
i  <  /  <  »  ♦  i 
l<4<A+ 1 

+  '*(06  —  (i6) 

fcfu*n  M< (input  value),  we  obtain  coordinate  positions  and  velocity  distribution 
m 

of  streamlines  by  using  the  calculation  results  of  the  printout.  Otherwise,  utiliz 
the  loosening  factor  r]  (input)  of  the  streamline  coordinates  to  obtain  new  stream¬ 
line  coordinates: 

f,  +  tj(re  —  rf)=»  'fl 

»,  +  »?(*  c  —  (17) 

01  +  t](0c  — 
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Again,  return  to  stop  5 


for  iterated  calculations  until  M<£  . 

in 

Now,  we  make  several  explanatory  notes  about  the  solution  process: 

1.  Just  as  for  procedures  in  report  (6),  when  solving  for  a,  first  rotate 

the  r-z  coordinates  by  45°  and  then  apply  the  method  of  obtaining  first  order  deriv¬ 
ative  by  using  sampling  curves. 

2.  Just  as  for  procedures  used  in  report  (12),  when  solving  for  r  ,  use  two 
sampling  curves  to  improve  the  precision  with  the  same  conditions  of  terminal 
points  as  in  report  (6) . 


3.  During  trial  claculation  used  in  the  paper,  for  the  time  being  the  varia¬ 
tions  of  hj^  and  ,\  are  left  aside.  This  is  feasible  in  compressor  calculations  of 
superchargers.  Variations  of  h!  and  A  will  not  cause  any  other  difficulties. 


4.  As  for  loss  revisions,  we  may  use  the  correction  method  of  total  pressure 
loss  in  report  (6) .  However,  for  convenience  this  program  adopts  a  correction  to 


the  multivariant  index: 


Here,  k  is  the  adiabatic  index,  n  is  the  multivariant  index;  and  r;  is  the  multi- 

o 

variant  efficiency.  The  value  of  nQ  is  determined  by  input. 


III.  Analysis  and  Discussion  of  Results 


A  program  was  compiled  by  using  this  method;  trial  calculations  were  conducted 
on  a  centrifugal  compressor  of  a  Number  780  supercharger  designed  by  the  Shanghai 
Shipyard,  the  Shanghai  Institute  of  Ship  Transportation  and  other  units.  Refer  to 
Figure  5  for  outlines  of  the  meridian  planes  of  a  compressor  blade  rotor.  The 
following  are  the  performance  parameters  thus  calculated:  inlet  temperature  T.=303°K; 

3  1 

inlet  specific  gravity  =  1.149  kg/m  ;  isopiestic  specific  heat  0^=1003.56  kg-m/kg; 
multivariant  index  n  =1.509;  flow  weight  kg;  angular  velocity 

w*880  radians/sec;  number  of  blades  Zn=20;  number  of  ^uasi  perpendi¬ 
cular  intersecting  planes  m=12;  number  of  flow  surfaces  (n+l)=9; 
number  of  S flow  surfaces  (h+l)=9. 

The  common  parabola  forming  method  is  used  to  generate  outlines  of  air-guiding 
rotors;  the  working  rotor  has  diametrical  straight  blades. 
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0  0.1  0.2'  *'w  " 

00  00 
Fig.  5 

(a)  Shapes  of  Meridian  Streamlines  and  Meridian  Plane  Outlines  of  a 
Blade  Rotor;  (b)  Distribution  of  Relative  Mach  Number  on  the  Averaged 
S2  Flow  Surface  (k=5)  . 

Key:  (1)  Meter  (m) ;  (2)  Quasi  perpendicular  intersecting  lines. 


(a)  Streamline  Shapes  on  a  Rotor  (j=l);  (b)  Distribution  of  Relative 
Mach  Number  on  a  Rotor  (j=l). 

Key:  (*)  Radian 


Fig.  9  Distribution  of  relative 
Mach  number  on  a  quasi  perpen¬ 
dicular  intersecting  line  at  the 
inlet . 

Key:  (1)  Rotor;  (2)  Stator. 

The  calculation  results  are  plotted  in  Figures  5-9.  Since  there  are  many  data, 
for  sake  of  simplicity,  these  curves  drawn  are  streamlines  of  the  trans-blade  flow 
field  for  j  =  1 ,  5  and  9;  the  corresponding  velocity  distributions  for  k=l,  3,  5,  7 
and  9;  the  streamline  distribution  for  k=5  (i.e.,  the  average  surface  of  flow 
channels);  and  the  velocity  distribution  for  j=l,  3,  5,  7  and  9. 

from  the  calculation  results,  we  drew  the  following  primary  conclusions: 

1.  We  can  see  from  Fig.  9  that  on  the  average  throughflow  surface 
(k=5),  the  distribution  of  relative  Mach  number  at  the  inlet  on 
the  quasi  perpendicular  intersecting  lines  is  quite  similar  to  that 
shown  in  Fig.  7  of  report  (12).  The  shapes  of  velocity  distribution 
and  streamline  plot  spectra  in  Figs  5-8  are  also  consistent  with 
curves  in  other  reports. 

-4 

2.  For  this  example,  whenfG  =  l°o,  £m=10  and  the  loosening  factor  r\  is 
selected  as  0.1,  the  three-dimensional  flow  field  solution  can  be  obtained  by 
applying  more  than  20  successive  substitutions.  Besides  the  printout  time,  the 
processing  time  is  approximately  8-10  minutes  by  using  our  university's  719  computer 
(125,000  operations  per  second).  This  processing  time  is  considerably  less  than 
the  mutually  successive  substitutions  of  nine  flow  surfaces  and  nine  S2  flow 
surfaces  in  several  calculations. 

3.  The  shortcoming  of  this  method  is  the  same  as  the  common  streamline 
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curvature  method:  it  is  generally  difficult  to  calculate  precisely  the  flows  at 
inlet  and  exit  --  corrections  are  required.  However,  the  revision  of  the  three- 
dimensional  flow  field  differs  from  the  situation  of  solving  separately  the 
and  S,  flow  surfaces.  The  work  will  be  continued  later  on. 

4.  Since  the  method  in  the  paper  is  based  on  the  assumption  that  Sj  flow 
surface  is  a  surface  of  revolution  with  more  limitations,  one  still  cannot  calcu¬ 
late  the  curving  state  of  the  flow  surface.  As  a  continuation  of  the  paper, 
later  this  assumption  will  be  eliminated  in  order  to  derive  a  real  solution  of  the 
three-dimensional  flow  field. 


APPENDIX 


In  the  following,  we  deduce  Equations  (4),  (6)  and  (7). 


Since  the  quasi  perpendicular  intersecting  plane  we  select  is  a  plane  of 
revolution,  the  generating  line  is  a  straight  line  (AB  in  Fig.  2a)  of  a  meridian 
plane: 


r  —  r,  +  qm  caxfi  i 
M  "  *,  —  qm  tin 


(Al) 


In  the  equation,  (r^,  z^)  are  coordinates  of  intersecting  point  c  between  the  blade 
rotor  and  the  quasi  perpendicular  intersecting  line  AB;  qm  is  the  intercept  along 
line  AB  starting  from  point  c;  and  \p  is  the  included  angle  between  z  axis  and  n^, 
the  normal  line  of  AB. 


By  eliminating  qm  in  Equation  (Al),  the  straight  line  equation  is 

(c  —  r.)  »in  ^  coi  •»  0  (A2) 


If  the  equation  of  blade  pressure  surface  is 

-  U'f*)  (A3) 

the  conditions  of  the  intersecting  curve  (space  curve)  equation  between  the  pressure 
surface  and  the  quasi  perpendicular  intersecting  plane  are 


u  •*  +  «.  «•> 

O  “  Iti'it  *<)  “  $<(*») 


(A4) 


The  direction  of  q  in  Equation  (1)  takes  on  the  direction  of  intersecting  curve. 


then 
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From  Equation  (A5) ,  we  know 


_ 1 '.i.  -  to>  4 

llm 


(A5) 


(A6) 


By  substituting  Equations  (A6)  and  (2)  into  Equation  (A5)  and  by  simple  algebraic 
manipulations,  we  obtain  Equation  (4). 


Since  along  the  arc  sector  of  qm=constant,  there  are 


it 

19 


o. 


19 

19 


I 


CA7) 


Substitute  q  in  Equation  (1)  by  9,  then  Equation  (6)  can  be  obtained  by  using 
Equation  (A7) . 


The  general  equation  of  flow  weight  streaming  into  any  curved  surface  S  is 


HIS 


If  the  equation  of  curved  surface  S  is  "  0 

then  the  unit  normal  vector  of  curved  surface  is 


VJ 

lvs| 


OS  . 


OS  OS  . 

roe  '*  +  a*  '• 


\vs\ 


(A8) 

(A9) 

(A10) 


In  the  present  situation,  the  curved  surface  S  is  quasi  perpendicular  inter¬ 
secting  plane,  the  equation  of  which  is  Equation  (A2) .  Therefore  its  unit  normal 
vector  is 


ii  ■  tin  -r  ji$  +  co»  4> 


(All) 


The  vector  form  of  relative  velocity  W  is 

W  -  H'.i.  i-  li  ,i,  +  H  .i, 


Hence, 


(A  12) 


W  •  l»  mt  If,  (in  <!>  +  H’,CO u/l 

™  If  a  tin  «  tin  4>  +  Wm  cot  m  cot  i/i 
—  If  a  cot  —  a) 

1 6 


(A  13) 


Besides,  the  area  micro -element  of  this  quasi  perpendicular  intersecting  plane  can 
be  written  as 


iS  -  r  J9  *4.  (All) 


Substituting  liquations  (A13)  and  (All)  into  liquation  (AS),  we  obtain  Equa¬ 
tion  (7). 


The  density  equations  arc  the  same  as  in  report  (7),  therefore  the  paper  does 
not  list  them. 


LITERATURE 


1.  Wu,  Zhonghua  A  General  Through-flows  Theory  of  Fluid  Flow  With  Subsonic  or 
Supersonic  Velocity  in  Turbomachines  of  Arbitrary  Hub  and  Casing  Shape,  NACA 
TN-2302  (1951). 

2.  Wu,  Zhonghua  A  General  Theory  of  Three-dimensional  Fluid  Flow  in  Axial  Flow, 
Diametrical  Flow  and  Mixed  Flow  Type  Subsonic  and  Supersonic  Velocity  Turbomachines, 
NACA  TN- 2604  (.1952). 

3.  Marsh,  H.  A  Digital  Computer  Program  for  the  Rectangular  Throughput  Flow 
Method  to  Calculate  Fluid  Flow  in  Arbitrary  Turbomachines,  ARC  R  5  M  3509  (1968). 

4.  Smith,  D.  T.  L.,  and  Frost,  D.  H.  Calculation  of  Fluid  Flowing  Past  Blades 
of  Turbomachines,  Proc.  IMF.,  Vol.  184,  Pt  3G  (I)  (1969-1970). 

5.  Bosman,  C.  Unsteady  State  and  Elimination  of  the  Need  for  Calculating  Fluid 
Flow  in  Turbomachines,  ARC  R  6  M  5746  (1974). 

6.  Katsanis,  T.  Use  Arbitrary  Perpendicular  Intersecting  Lines  to  Calculate 
the  Distribution  of  Fluid  Flow  on  Meridian  Plane  of  Turbomachines,  NASA  TN-D  2546 
(1964) . 

7.  Katsanis,  T.  Use  Arbitrary  Perpendicular  Intersecting  Lines  to  Calculate 
Distribution  of  Fluid  Flow  on  Trans-blade  Surface  of  Turbomachines,  NASA  TN-D  2809 
(1965) . 

8.  Katsanis,  T.  Use  Arbitrary  Perpendicular  Intersecting  Lines  to  Calculate  the 
Distribution  of  Fluid  Row  in  Turbomachincs,  AJME,  Ser.  A,  88,  2  (1966). 

9.  Novak,  R.  A.  A  Radius  of  Curvature  Calculation  Method  for  Solving  Fluid  Flow 
Problems,  ASME,  Ser.  A,  89,  4  (1967). 

10.  Meiwei,  Taili  [Chinese  transliteration  of  characters  for  Japanese  name]  and 
Zhonglai  Jingchi  [Chinese  transliteration  of  Japanese  name]  Rotor  Blade  Theory  of 
Arbitrary  Surface  of  Revolution,  ASME,  Ser.  A,  93,  4  (1971). 

11.  Katsanis,  T.  A  Fortran  Program  for  Velocity  Distribution  and  Blocking  Flow 
at  Surface  of  Turbine  Blade  Array,  NASA  TN-D  6177  (1971). 

17 


END 

* 

DATE 

FILMED 


DTIC 


